PROCEEDINGS 

OF SCIENCE 



o 

(N 
-i— > 
O 

O 



Oh! 



Gluon Mass in Landau Gauge QCD 



P. Bicudo* 

CFTP, Departamento de Fisica, Institute Superior Tecnico, Av. Rovisco Pais, 1049-001 Lisboa, 
Portugal 



E-mail: bicudo@ist . utl ■ pt 



0. Oliveira 

Departamento de Fisica, Universidade de Coimbra, 3004-516 Coimbra, Portugal 



E-mail: orlando@teor. fis.uc.pt 



> 

in 

On 

d 
o 



The interpretation of the Landau gauge lattice gluon propagator as a massive type bosonic propa- 
gator is investigated for i) an infrared constant gluon mass; ii) an ultraviolet constant gluon mass; 
iii) a momentum dependent mass. We find that the infrared data can be associated with a mas- 
sive propagator with a constant gluon mass of 651(12) MeV, but the ultraviolet lattice data is not 
compatible this type of propagator. The scenario of a momentum dependent gluon mass gives a 
decreasing mass with the momentum, starting from a value of ~ 630 MeV in the infrared region 
and suggesting a q 2 \nq 2 dependence for momenta above 1 GeV. 
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1. Introduction and Motivation 

If at the classical level SU(3) Yang-Mills theory is conformal invariant, the corresponding 
quantum theory breaks this symmetry via dimensional transmutation. Indeed, dimensional trans- 
mutation introduces a scale in QCD, Aqcd ~ 300 MeV, which defines the typical energy for strong 
interactions. Despite the breaking of conformal invariance, at the level of the lagrangian a mass 
term is forbidden by gauge invariance. Moreover, within the perturbative solution of QCD in Lan- 
dau gauge the gluon is a massless particle. However, if one goes beyond perturbation theory and 
looks for nonperturbative solutions of theory, then a dynamical mass for the gluon becomes possi- 
ble [j|]. This gluon mass is, typically, a function of the gluon momentum M(q 2 ). 

A non-vanishing gluon mass is welcome to regularize infrared divergences and solve some 
problems related with unitarity. Diffractive phenomena ^ and inclusive radiative decays of J/yr 
and T [J3j] suggest a massive gluon with a mass in the range 0.500 - 1.2 GeV depending on how you 
define the mass. Moreover, lattice simulations also suggest an infrared gluon hard mass of ~ 600 
MeV [|] and an ultraviolet mass M g ~ 1.0 GeV [g, g]. 

A massive gluon is also welcome within the dual picture of the QCD vaccuum, where a Meiss- 
ner effect due to chromomagnetic Abrikosov flux tubes introduces an effective gluon mass [J7|, [| Pj. 

From the point of view of the Dyson-Schwinger equations, the idea of a gluon mass which as 



a function of the gluon momenta fits naturally within the so-called decoupling solution [10, 11, 12, 
13, 14]. Indeed, the numerical solutions of the DSE give a M(q 2 ) which takes its largest value at 



zero momentum, where M (0) ~ 600 MeV, and vanishes for q S> Aqcd, recovering, in this way, the 
usual perturbative propagator at high momentum. 

Lattice QCD simulations also provide support for a non-vanishing gluon mass - see, for ex- 
ample, [g, g, |, 15] and references there in. The precise value for M(q 2 ) depends on how the gluon 



propagator is modeled. 

2. The Gluon Propagator and The Gluon Mass 

In the Landau gauge the momentum space gluon propagator is given by 

D%(q 2 ) = (A^q)A b v (-q)) = 8 ub (V - ^) D(q 2 ) . (2.1) 
As definition of the gluon mass we take 

and, in the following, we will measure Z(q 2 ) and M 2 (q 2 ) from the gluon data computed in SU(3) 
lattice QCD simulations at j8 = 6.0 for various volumes. Details of the simulation can be found in 



[|15J]. 

The gluon propagator used to compute M 2 {q 2 ) are renormalized at /I = 3 GeV such that 
D(jU 2 ) = 1/jU 2 - see [JT3J] for details. 

The lattice spacing effects are removed/reduced applying, for each volume, a conic cut for 
momenta q > 1 GeV [S]. To achieve a good description of the infrared region, for momenta smaller 
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Renormalized Gluon Propagator - Jl = 3 GeV 

P = 6.0 - combined 32*, 4S 4 , 64 4 , 80 4 data 
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Renormalized Gluon Dressing Function - Jl = 3 GeV 

P = 6.0- combined 32 4 , 48 4 , 64 4 , 80 4 data 
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Figure 1: Renormalized gluon propagator D(q ) (left) and gluon dressing function q D(q ) (right). 

than 1 GeV all the data points coming from the various simulations were considered. Further cuts 
were applied to reduce finite volume effects. The data from the smaller lattices which deviate from 
our largest volume simulation ~ (8.1 fm) 4 were deleted from the data set. Then, from the 64 4 
propagator only data with q > 425 MeV was considered, from 48 only q > 671 MeV data was 
considered and from 32 only data with q > 848 MeV was included. However, the IR and UV 
cuts described are not able to remove all lattice spacing and volume effects. For each lattice and 
for the same q 2 coming from different q^, if the different estimates of the propagator don't agree 
within one standard deviation, one of the points is removed. For example, for the 80 4 lattice for 
momentum q = 457 MeV there are two estimates for the gluon propagator, D(q 2 ) = 6.515(96) 
GeV 2 and D(q 2 ) = 6.382(64) GeV -2 , coming from different types of momenta. The first value 
is clearly above all the data points and it was not considered in the combined data set. In this way, 
the surviving points will produce a unique curve for D{q 2 ). 

The combined lattice data for gluon propagator, after all the cuts, is shown in figure [T] together 
with the dressing function q 2 D(q 2 ). 

3. An Infrared Constant Mass 

Let us start our discussion on the gluon mass considering the case of an infrared constant gluon 
mass, i.e. assuming that the gluon propagator is described by 



D(q 2 



q 2 +M 2 



(3.1) 



where Z and M are constants, in the momentum range [0,q ma x]- The results of fitting the lattice 
data to (3. 1 ) are plotted in figure |[ 

Figure || shows a M 2 and Z that are, within one standard deviation, stable against variation of 
qmax- Furthermore, demanding that % 2 jd.o.f . < 1.8, then the infrared propagator can be described 
by equation ( 3. 1 ) for momenta up to q ~ 530 MeV. The fits reported here have a 1 .26 < % 2 / d.o.f. < 
1.8. 

From the above results one can conclude that, with the possible exception of the deep infrared 
region, D(q 2 ) can be described by a massive type propagator with a constant mass in the low energy 
regime. 
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Figure 2: Fittins an infrared constant mass (Y 2 /d.o.f. < 1.8). 




Figure 3: Lattice gluon propagator and the infrared fit the smallest % 2 /d.o.f. 

The fit with the lower X 2 /d.o.f. has q max = 466 MeV, Z = 4.05(10) and M = 651(12) MeV. 
The fit together with the lattice propagator are reported in figure |3[ 

4. An Ultraviolet Constant Mass 



If one applies the same reasoning to the ultraviolet (UV) region, i.e. for q > 2.5 GeV, it 
turns out that M 2 is not stable against variation of the range of momenta considered. In this sense 
one cannot define a gluon constant mass to the ultraviolet region. This is not a surprise. Indeed 
to renormalize the lattice propagator one uses a perturbative inspired one-loop expression, which 



describes very well the data in the UV region - see []15j] for details. 

The results discussed here for the UV are not in contradiction with thise of [Q, |6|, where an 
ultraviolet gluon mass of ~ 1 GeV was claimed. In [Q, ||] an ultraviolet regulator was used and the 
full set of lattice data surviving the conic cut fitted to 



D(q 2 ) 



\\n(q 2 +M 2 )(q- 2 +M- 2 ) 



q 2 +M 2 



(4.1) 



where M is the gluon mass. The difference is due to different definitions for the gluon mass. 

In conclusion, in what concerns the ultraviolet region, the lattice gluon propagator is not de- 
scribed by a constant massive propagator. 



4 



Gluon Mass in Landau Gauge QCD 



P. Bicudo 



Renormalized Inverse Gluon Propagator - (i - 3 GeV 

P = 6.0 - combined 32*, 48 4 , 64*. 80 + data 




Figure 4: Inverse gluon propagator - note the "almost" linear behaviour with q 2 of 1 /D(q 2 ). 
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Figure 5: Overall picture for running gluon mass - only data whose fit has x 2 /d.o.f. < 1.8 and relative error 
less than 30% is included in the plot. 

5. A Momentum Dependent Massive Gluon 



In this section we will assume a D{q 2 ) given by equation (2.2). The momentum functions 
M(q 2 ) and Z{q 2 ) will be referred as the running mass and running dressing function, respectively. 

For the computation of Z(q 2 ) and M 2 (q 2 ) two different methods called (i) linear 4 points 
and (ii) massive 4 points in the figures, were considered. Furthermore, Z(q 2 ) and M 2 {q 2 ) were 
investigated using both D{q 2 ), see figure [j], and \/D{q 2 ), see figure |j. The interested reader can 
find the details on the procedure to extract the functions from the lattice data in [O]. The fits with 
a % /d.o.f. > 1.8 were exclude from the analysis. 

The running mass, computed with both methods, is shown in figure ||. M 2 (q 2 ) is positive 
in the infrared region, decreases with q and becomes negative around q ~ 800 MeV. Although 
M 2 (q 2 ) becomes negative, q 2 +M 2 (q 2 ) is always positive defined, i.e. the propagator has no poles 
for euclidean momenta. 

For the infrared region, the running mass measured for the smallest momenta have q = 249 
MeV (method i) and 273 MeV (method ii). The corresponding mass values are, respectively, 
652(79) MeV and 669(88) MeV. We note the excellent agreement with the estimation of a hard 
infrared constant mass 651(12) MeV (see section ||). 

Let us investigate a possible ansatz for M 2 {q 2 ). Our best fit occurs when the ultraviolet region 
and the infrared region are studied separately. Given that the statistical errors on M 2 (q 2 ) increase 
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Figure 6: Running gluon mass M (q ) and the fits (full lines) to (5.1) and (5.2). 



with q, it compromises the investigation of the ultraviolet behavior - see figure ||. However, starting 
at a relatively low momenta, let us say around 1 GeV, one can test for the q 2 dependence of M 2 (q 2 ). 
Our best fit points towards a M 2 (q 2 ) ~ q 2 lnq 2 at high momentum and an infrared dependence 
where M 2 (q 2 ) ~ q 2 . 

The outcome of the separate fits to the two momentum regions is 



M 2 (q 2 ) = 



0.534(15) - 0.943(48) g 2 
0.578(37)- 1.239(64) q 2 



method (i) 
method (ii) 



(5.1) 



with a % 2 /d.of. of 0.6 and 1.6, respectively, for the infrared region, i.e for q < 1 GeV, and 



-0.349(22) -0.1465(79) ^ 2 ln^ 2 , method (i) 
-0.288(33) - 0.2050(78) q 2 Inq 2 , method (ii) 



(5.2) 



with a % 2 jd.of . of 0.5 and 1.8, respectively, for the UV region, i.e. for q > 1 GeV. In the above 
formula M 2 (q 2 ) and q 2 are given in GeV 2 . 



In figure owe plot the lattice data together with the fits to (5.1) and (5.2). The fits to the infrared 



region, see equation (pT|), give an M(0) 731(11) MeV (method i)and 760(24) MeV (method ii), 
which is slightly larger than the constant infrared mass computed in section ||} Anyway, the fits 
to the infrared region suggest a finite D(0). Moreover, the predicted M(0) is associated with a 
D(0) ~ 2 GeV" 2 , a value of the same order of magnitude as predicted by recent large volume 



lattice simulations [16] 



The running gluon dressing function Z(q 2 ) data is displayed in figure ^ together with the dress- 
ing function the (full line) computed using the 1-loop perturbative expression for Z{q 2 ) rescaled by 
0.7. 

In the low momentum region, Z(q 2 ) decreases from q = up to 1 GeV. However, when the zero 
momentum is approached from above, Z(q 2 ) seems to saturate around q = 400 MeV. Unfortunately, 
the large statistical errors in Z{q 2 ) for q < 1 GeV make it difficult to disentangle the functional 
dependence but, clearly, Z(q 2 ) does not follows the perturbative behavior. 

The gluon dressing function Z(q 2 ) is well described by the anstaz 



2 ) 



[A + ln(<7 2 + m 2 )] 



(5.3) 
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Figure 7: Running gluon dressing function Z(q 2 ). Only data whose fit has a x 1 jd.o.f . < 1.8 and relative 
error less than 30% is included in the plot. For the meaning of the full line, see the main text. 

where y = 13/22 is the anomalous gluon dimension, over the full range of momenta. The fits 
give Z = 1.048(75), A = -0.43(21), ra 2 = 1.57(33) GeV 2 for a % 2 /d.o.f. = 1.8 (method i) and 
Z = 0.98(14), A = -0.54(38), = 1.76(68) GeV 2 with a % 2 /d.o.f. = 1.9 (method ii). 
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